in 

o 
o 

(N 



^ 



X 



The pole structure of the unitary, crossing 

symmetric low energy tttt scattering 

amplitudes 

Z. Y. Zhou, G. Y. Qin\ P. Zhang, Z. G. Xiao and H. Q. Zheng 
Department of Physics, Peking University, Beijing 100871, P. R. China 

and 

N. Wu 
Institute for High Energy Physics, Chinese Academy of Science, Beijing 



^. : 100039, P. R. Chma 

(3 ■ February 2, 2008 

(N' 

Abstract 

(^ 

^ I The pole structure of the low energy vrvr scattering amplitudes is 

studied using a proper chiral unitarization method combined with cross- 

pg . ing symmetry and the low energy phase shift data. It is found that the 

\0 ; a pole position is at M^ = 470 ± 50MeV, T„ = 570 ± bOMeV. The 

^-p ' existence of the virtual state pole in the IJ=20 channel is reconfirmed. 

^^ . Various threshold parameters are estimated and are found in general 

"^ I in good agreement with the results obtained from the Roy equation 

p i' analyses. 

D , Key words: tttt scattering; unitarity; dispersion relations, a resonance 



PACS number: 14.40.Ev, 13.85.Dz, 11.55.Bq, ll.30.Rd 



The a or the /o(600) resonance has attracted much interest in recent 
years P] and more and more physicists, if not all, believe that it truly exists as 
a light and broad resonance. [2] There were two main difficulties to accept the 
a resonance, one is that being light and broad makes the a's contribution to 
the phase shift hard to be distinguishable from background contributions while 
the latter are often out of control in various phenomenological studies. How- 
ever, a careful analysis revealed that the cut contribution to the phase shift 
is negative in the I,J=0,0 channel tttt scattering [3], and thus the a resonance 
has to exist to saturate the experimental data. Another objection is that the 
TTTT phase shift in the 1=0, J=0 channel does not pass 7r/2 at moderately low 
energies which should have appeared if there exist a light a as required by the 
standard Breit-Wigner description of resonances. However it is pointed out 
that a light and broad resonance rejects a simple Breit-Wigner description. As 

"'^Present address: Department of Physics, McGill University. 3600 University Street, 
Montreal, Quebec, H3A 2T8, Canada. 



a consequence a light and broad resonance's contribution to the phase shift 
may pass 90° only at energies well above IGeV. jl] 

Although the existence of the a resonance seems to be firmly established, a 
precise determination to the location of this resonance on the complex s plane 
can still be challenging. A summary of the previous determinations on the pole 
location of the a resonance may be found in Ref. [H] and the results are found 
to be rather diverse, which is reflected in the Review of Particle Properties as 
the mass of the a ranges from 400MeV to 1200MeV. In more recent studies, the 
pole mass of the a resonance is estimated to be M ~ SOOMeV and F ~ 640MeV 
in Ref. [Sj, in Ref. [H] an extensive analysis using Roy equations indicates 
M = 470±30MeV, F = 590±40MeV. On the other side, the E791 experimental 
group claim from the D^ -^ Tr+Tr^TT"*" data analysis that M = 478^23 =tl7MeV, 
F = 324l;4Q ± 21MeV,^ [7j whereas the BES collaboration has found from 
BES-II J/ijj -^ unir data the a pole mass to be M = 541 ± 20 ± 32MeV, 
F = 504 ± 50 ± 60MeV jS] and BES-I J/tp -^ umr data the pole mass as 
M = 434±78MeV, F = 404±86MeV. Under this situation it is worthwhile 
to make further efforts towards a more accurate determination to the a pole. 
What we will adopt here is a new approach recently proposed by our group 
which is elaborated in details in Ref. ^U] discussing ttK scatterings. The 
approach respects all known fundamental properties of S matrix theory except 
that crossing symmetry is not implemented automatically. We will remedy 
the deficiency here by considering the constraints from crossing symmetry, 
the so called Balanchandran-Nuyts-Roskies (BNR) relations in the fit. In the 
following we begin by briefly review the new unitarization scheme developed 
in Ref. [10]. 

For a partial wave scattering in a given channel, the physical S matrix can 
in general be factorized as, 

^Pfty _ -Q ^ft . ^c«t ^ (1) 

i 

where S*^* are the simplest 5* matrices characterizing the isolated singularities 
oiSP^y ^, which are: 



1. for a virtual state pole at s = Sq (0 < Sq < 4m: 



■K>1 



1 — tp[s)a 

The scattering length contributed by a virtual state pole defined as such 
is 

(3) 



4m2 - So 



^Corresponding to a pole mass: M = 489MeV, T = 346MeV. 



2. for a resonance located at zq (and Zq) on the second sheet, the S matrix 
can be written as, 

M^\zo] - s + ip{s)s /"^["°' 



M2[zo] -s-zp(s)s 



Im[zo] 



where 



lm[^yzo{zo -4m: 



Re[^zo{zo - 4m2)] 



(5) 



The expression of the scattering length contributed by a resonance defined 
as above is also easily obtainable. The S''^"* contains only cuts which can be 
parameterized in the following simple form, 

TT J R [S' - Sq)[s' - S) 

where L = (—00, 0] and R denotes cuts at higher energies other than the 27r 
elastic cut. To a good approximation R starts at 4m\-. The parameterization 
form, Eqs. (Q) - ©) which may be viewed as an improved version of the 
Dalitz-Tuan representation, has the advantage that S"^' does not contribute 
to the discontinuity of /,Pn] which means the following relation: 

Im^,^/(.) = -^log|S^^n^)|. (8) 

The left cut integral in the above equation are estimated by replacing S*^^^ 
on the r.h.s. of Eq. (jHI) by S^^'^ which is available at two-loop order for tttt 
scatterings ^5. In Fig. Qwe plot the function ImL/(s) using both tree, 1-loop 
and 2-loop xPT results at moderately low energies from which one may get 
an impression on the convergence speed of the chiral expansion. Evidences 
exist in the I,J=2,0 channel of vrvr scatterings [3] and I,J=|,0 channel of nK 
scatterings jTU] to support such an approximation. The present method is cer- 
tainly superior to the method of Ref . j2] , since now the constraint of unitarity 
is automatically embedded in the formalism. As a bonus of being unitary, 
when evaluating the left hand integral using the following equation, 

ImL/(s) = -^log|5^^^(.)|, (9) 

the once-subtracted dispersion integral automatically converges owing to the 
logarithmic form of the integrand. This property further regulates and sup- 
presses the bad high energy behavior of the perturbation amplitudes. Never- 
theless it should be pointed out that though the dispersion integral is formally 
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Figure 1: Tree, 1-loop and 2-loop results of Im2,/(s) from Eq. 0: left) 1=0, 
J=0; right) 1=2, J=0. The low energy parameters of the O(p^) chiral La- 
grangian are taken from Ref. 



convergent when one uses xPT result to make numerical estimation, the con- 
tribution at high energies is beyond control theoretically due to the bad high 
energy behavior of xPT. Therefore we introduce a cutoff parameter when eval- 
uating the left hand integral in each channel in the numerical studies, which are 
determined by the minimization of x^- This is of course, only a phenomenolog- 
ical approach. Actually, the above parametrization form can be obtained only 
by assuming the analyticity property on the whole cut plane, which is a conse- 
quence of Mandelstam representation. The Mandelstam representation how- 
ever goes beyond what is rigorously established from field theory. Nevertheless 
the analyticity domain is large enough for phenomenological applications not 
very far from the physical region. ^21 

The above discussion briefly reviewed what is already established in previ- 
ous studies. However further improvement is possible. In tttt scatterings the 
kinematics is much simpler than in -kK scatterings which makes it easy to re- 
alize that the subtraction constant appeared in Eq. (Q) can actually be fixed. 
To see this we recast Eq. (^ into the following form: 

^log(l + 2zp(.)T^^'n.)) = ^5^1og(^^'(^)) + /(^) • (10) 

The l.h.s. of the above equation contains the normal branch points of 5* as 
well as the branch points at Zi (which denote the pole positions of 5). On the 
r.h.s of Eq. ()1U|) the normal branch points are presented in /(s) whereas the 
branch points at Zi are contained in the first term. In the absence of the long 
range Coulomb interaction the vrvr scattering amplitudes T"^^^ is well behaved 
when s — > and 4m^. Also it is not difficult to determine the behavior of 
S^^{s) when s — > and 4m^. As a consequence, one obtains, by evaluating 
the asymptotic behavior of each side of Eq. (fTUj) when s — > and s — >• 4m^, 
the following relations: 



/(O) = 
a = y ^a 



/(4m^) 



(11) 
(12) 



respectively. In Eq. (I12|l a is the scattering length of the partial wave ampli- 
tude, a^' denotes the scattering length contributed by pole Ri. One advantage 
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Figure 2: Different estimates on the baclcground contribution to the phase 
shift. Left: in the 1=2, J=0 channel; right: in the 1=0, J=0 channel (here the 
right hand cut contribution is not included). The experimental data of 6q are 
from Ref. 



of using Eq. (fTTj) is that it fixes the subtraction constant appeared in Eq. ((7j) 
and suggests that the natural subtraction point in Eq. (jTj) is at Sq = 0, that is. 



f{s) 



.^ s'is' - s) 



■ds' 



TX 



i s'{s' - s) 



ds' 



(13) 



The validity of Eq. (fT^ can also be demonstrated by directly evaluating the 
asymptotic behavior of f{s) when s -^ 0. Starting from Eq. p3|) and Eq. (jH)), 
remembering that T^'^s^(O) is a constant, it is not difficult to demonstrate that 
the asymptotic behavior of S''^"*(s) when s — > agrees with the l.h.s. of Eq. (P).^ 
Therefore Eq. (jl3|) is correct. 



The Eq. (^ has been found very useful in Ref. fU] when studying the 
low lying resonance structure in ttK scatterings. Motivated by the discussion 
given earlier in this Letter it is worthwhile to re-examine the a pole in tttt 
scatterings. For such a propose, it is necessary to firstly estimate the the 
background contribution /. The left hand cut integral in Eq. (J13p is estimated 
using two loop chiral perturbation theory results on Iuil/. The integral interval 
of the left hand integral in Eq. ()13l ) is taken as [— A|, 0]. In Fig. |21we plotted 
the cutoff dependence of the function /. For the right hand cut integral in 
Eq. (|13p. we neglect its contribution in the IJ=11 and IJ=20 channels. In the 
IJ=00 channel we will use a simple model found in the literature, ^H] which 
gives a rather good fit to the inelasticity parameter in the energy region above 
KK threshold, to estimate its contribution. The integration will be truncated 
at some appropriately chosen cutoff scale: A^ ~ 2GeV^. 

In the following we will discuss the speciality of IJ=20, 11, 00 channels 
separately. Then we will include the constraints from crossing symmetry to 
make a combined fit to the data of the three channels below IGeV and the so 
called BNR relations. 

IJ=20 channel: 



^Up to a sign if the total number of bound or virtual states in a given channel is odd. 



According to Eq. ()12|). if there is no resonance or virtual state in the IJ=20 
channel the scattering length will have a simple relation: a^ = /(4?7i^), which 
is calculable and its value depends only on A^ in our scheme. However such a 
formula gives a too large magnitude of the scattering length Qq. For example, 
when A| changes from O.SGeV^ to lO^GeV^, /(4m^) varies from —0.107 to 
—0.187, a value too large comparing with the experimental value Oq = — 0.028± 
0.012 ^n] and the lattice results ^Tj, though the lattice method is still not 
mature yet. The apparent discrepancy can be solved, recalling that there 
actually exists a virtual state pole in the IJ=20 channel. ^H] The pole location 
can be estimated from the 2 loop perturbative amplitude from xPT to be 
So = 0.049m^ on the real axis on the second sheet of the complex s-plane. 
According to the definition of pole term adopted in the present paper, this 
virtual state pole's contribution to Oq is rather large (~ 0.11) which can not 
be neglected. The existence of the virtual state pole is unavoidable if chiral 
perturbation expansion works in the small s region (actually it is a prediction 
of current algebra and the relativistic kinematics ^Hl)- In the present approach 
its existence is further tested and is found to be fully consistent with the S 
matrix unitarization of xPT. Further, the existence of such a state raises an 
interesting question on the deeper understanding of the physical interpretation 
of 5* matrix poles and the quark hadron duality. 

IJ=11 channel: 

The speciality in this channel is that the first term of the scattering length 
expanded in powers of pion momentum vanishes: 

1 



ReTuis) = q''[a'j + b^ + 0{q')] , {q=-^s-4ml), (14) 

which leads to the following constraint: 

f{^ml) = -J2a^^. (15) 

i 

It is however not convenient to solve this constraint to reduce one resonance 
parameter, we therefore make the trick to subtract / both at s = and at 

s = 4m^, 

f(s) = li^a^s . sis -Ami) f Im/(s') 



Ami TT Jl+rs'{s' -Aml){s' - s) 



E 



^R.^ sjs-Aml) r Imfjs') 

^rnl vr Jl-,r s'{s' - Aml){s' - sf' '^ ^ 



to meet both the threshold and the s = behavior of f{s). ^ There exists 
apparent disagreement between the phase shift data [ini EOj and theory in 
this channel, which results in a large total xj=ij=i- The contribution to the 

'*We are in debt to I. Caprini for pointing out the mistake we made here in the previous 
edition. 



large xli comes mainly from low energy region {^/s ~ 0.55GeV) and also from 
high energy region {^/s ~ O.QGeV), considering the extraordinarily small error 
bars of the data. As a consequence, the Xii is very sensitive to the physical 
fit parameters as well as A| in this channel. In addition to the p(770) pole, 
if we introduce another pole p' to simulate the effects of p(1450) and higher 
resonances it will be helpful to reduce the large total Xn by roughly 50 and 
the x^]^ is no longer very sensitive to A|. Also it will improve the fit value 
of the scattering length and the effective range parameter defined in Eq. fll4|) . 
Therefore we will add a p' pole contribution in the fit. ^ The resulting Xn is still 
very large mainly due to the apparent disagreement between the theoretical 
curve and the data at low energies. Nevertheless the p pole position is very 
stable against different choices of A|. Also the a pole position changes very 
little with or without adding the p' pole (less than 10 MeV) in the combined 
fit as described later. 

Because the consistency problem between data from Ref. ^31201 and theory 
as mentioned above, it is necessary to look for experimental information from 
other sources. There exists accurate low energy data on pion form-factor in 
the 1=1, J=l channel from the CMD-2 experiment. j23^ However since what 
we need is the phase shift data, the data provided by CMD-2 experimental 
group are not directly usable here. The way we extract the 'experimental' 
data on 6l is the following: we use Eq. (8) in Ref. [2Ij to refit the form-factor 
data to determine parameters, then we turn off the p-u mixing parameter 6 
and calculate the phase shift 6l using again Eq. (8) of Ref. j23- ^^ this way 
we can manipulate the 'experimental' data on 5[. It is arbitrary to choose 
the center of mass energy of each data point, nevertheless we fix the center of 
mass energies of the 'data' set the same as the original data for form- factors 
provided by Ref. [21]. The error bars of the phase shift 'data' are calculated 
using the error matrix method. In this way we obtain 43 'data' points for (5|. 

If using our formula to fit the manipulated 6l data described above we get 
a very small total Xn? which suggests that there may exist strong correlations 
between the error bars of the manipulated data and/or that our formula and 
the Eq. (8) of Ref. [21] are highly consistent. Though there appears drastic 
difference on total x^ between the data of Ref. JH] and the data manipulated 
from Ref. [2^ , it will be shown later that the difference has only small influence 
on the determination of the a pole position. 

IJ=00 channel: 



^When introducing the additional pole one has to confine it in the vicinity of 1450MeV. 
Otherwise the pole will move down to M ^ 300MeV with a large width to remedy the appar- 
ent disagreement between theoretical curve and data at low energies. This is unacceptable 
because the scattering length and the slope parameter will be totally different from the xPT 
results. The p(1450) pole can not be fixed reasonably well in the present calculation since 
it suffers from the fact that we do not attempt to include the data around 1450 MeV, the 
deficiency of the data at low energies, and uncertainty of the cut-off integral when evaluating 
the background contribution. Therefore to introduce the p' pole only aims at improving the 
present fit results. 

^We thank Prof. Leutwyler for bringing Ref. [23 to our attention. 



Similar to the argument given in Ref. [3], the cut contribution calculated 
from Eq. ()13|) is negative, which clearly demonstrates that the a resonance 
must exist to reproduce the observed experimental phase shift in the 1=0, 
J=0 channel, which is positive and rising at low energies. The determination 
of the pole position of the a resonance is straightforward. In the 1=0, J=0 
channel we begin by introducing two poles below IGeV region: one is the 
a resonance under investigation and the another one corresponds to /o(980). 
With their masses and widths as fit parameters, MINUIT displays an obvious 
broad resonance and a narrow resonance. However, cautious readers may still 
ask the question whether the strong enhancement of 5q below IGeV is owing 
to the sum of residue effects of many resonances located far away from the 
physical region other than from a single a pole. To explore such a possibility 
we introduce a third pole without any prejudice on where it should locate. 
We will test whether the a pole location will be severely distorted by the 
introduction of this additional pole. The answer is no. The a pole moves only 
slightly with or without the third pole included (roughly less than 20 MeV). 
On the other side the location of the third pole is found to be very flexible 
when varying A^, but always stays far away on the left. Roughly we find M3 
ranges from a few tens of MeV to a few hundreds of MeV whereas Fs ranges 
from 1 - 3GeV. Such a pole would better be described as pure background and 
its contribution to the phase shift remains small. 

The combined fit with constraints from crossing symmetry 
A major shortcoming of the unitarized partial wave amplitudes as described 
above is that crossing symmetry is not automatically built in, despite of its 
apparent advantages. On the other side, necessary conditions involving only 
low partial waves can be derived from crossing symmetry. In tttt scatterings, 
there exists crossing relations between partial wave amplitudes - BNR rela- 
tions. There exists five such relations if only s-wave and p-wave are considered, 
i.e., 



I: 



I [s-A){?,s-A)[tl{s) + 2tl{s)]ds = Q, 
Jo 

II: /" (s-4)i?°[2t°(s)-5t2(s)]rfs = , 

Jo 

III : {s- 4)i?°[2t[!(s) - 5tl{s)]ds - 9 / (s - 4)2/?JtJ(s) = , 

Jo Jo 

pA pi 

IV : I {s- A)Rl[2tl{s) - htl{s)]ds + 6 / (s - AfR\t\{s) = , 

Jo Jo 

V : f {s- 4)i?0[2t0(s) - 5tl{s)]ds - 15 [ (s - AfRltHs) = ,(17) 

Jo Jo 
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Figure 3: Fit results on phase shifts from Eq. (J2H)- Data in the IJ=11 channel 
are: open circle from Ref. jIH], full square from Ref. [201 • Data in the other 
channels are described in the text. 



where R^' are, 



Rq — 1 , Rq — I 



i?o = 3s - 4 , i?} = 5s - 4 , 

i?o = lOs^ - 32s + 16 , Rl = 21s2 - 48s + 16 , 

Rl = 35s^ - ISOs^ + 240s - 64 . (18) 

Dimensional parameters in Eqs. (J17II18|) are expressed in unit of ?7i^ for sim- 
plicity. All these BNR relations can be rewritten in a compact form: 



I 



ds[PUs)T^{s) + Pl,is)Tlis) + PUs)T^{s)] = 



where i ranges from I to V and Pjjis) are polynomials. We will make a 
combined fit on the IJ=20,11,00 channels together with the 5 BNR relations 
using the penalty function method. |22] We define the penalty function as 
follows: 



i=i 



.2 2-^ 



J^'ds[PUs)TS{s) + Ph{s)Tl{s) + Piois)Ti{s) 



1 llo ds{\PMTo%s)\ + \Pl^{s)n{s)\ + \Pio{s)T,\sm^ 

(19) 

and define the total x^ as, 

Xtot = X20 + xli + Xoo + xInr ■ (20) 

Altogether we have 10 physical parameters corresponding to 5 poles (3 in the 
IJ=00 channel and 2 in the IJ=11 channel). In addition we have three cutoff 
parameters when evaluating the background contributions in three channels. 
The penalty factor, e~^, varies between 10~^ to 10~^. The rule to choose such a 
range is that we require that the violation of each BNR relation (corresponding 
to the value of Err^^^ defined in Eq. (I19|) ) is at the level of one percent. 

In the following combined fit we will test different combination of data. 
Firstly we use the data of IJ=11 channel from Ref. f^l 120] ■ Below 1 GeV 
there are 23 data points in the IJ=20 channel ^3], 39 data points in the IJ=11 
channel dH 120] and 36 data points [I1I221 123] in the IJ=00 channel. In the 
fit it is found that the total x^ is sensitive to (A|)/j=2o and the best choice 
corresponds to (A^)/j=2o — 3.4GeV^. On the other side the total x^ is no longer 
sensitive to (A|)/j=ii when the p' pole is included. Also the total x^ and all 
major physical outputs are insensitive to (A|)7j=oo- Part of the reason behind 
such a phenomenon is because in the IJ=00 channel the third pole's location 
is sensitive to the choice of the cutoff parameter and it absorbs the cutoff 
dependence of the background integral. It is found that the a pole location is 
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Our results 


xPTm 


Roy Eqs. 


Exp. pi 


Unit 


< 


0.211 ±0.011 
0.264 ±0.015 


0.220 ±0.005 
0.280 ±0.011 


0.220 ±0.005 
0.276 ±0.006 


0.26 ±0.05 
0.25 ±0.03 


m;2 


bl 


-0.440 ±0.011 
-0.785 ±0.010 


-0.423 ±0.010 
-0.762 ±0.021 


-0.444 ±0.010 
-0.803 ±0.012 


-0.28 ±0.12 
-0.82 ±0.08 


10-1 
lO-im-2 


b\ 


0.367 ±0.003 
0.563 ±0.003 


0.380 ±0.021 
0.58 ±0.12 


0.379 ±0.005 
0.567 ±0.013 


0.38 ±0.02 


lO-im-2 
10-2m;4 



Table 1: Threshold parameters comparing with other works. Phase shift data 
in the IJ=11 channel are from Ref. fHll2n]- More accurate experimental results 
of Oq and Og of most recent KeA experiment by E865 Collaboration can be found 
in reference |21I : a[j = 0.216 ± 0.013(stat.) ± 0.004(syst.) ± 0.002(theor.), 



al = -0.0454 ± 0.0031 (stat.) ± 0.0010(syst.) ± 0.0008(theor.) . See also the first 
reference of |2H] for a different analysis of the same data. 

almost inert against the variation of the cutoff parameter. Though the third 
pole's position is very flexible, it always wanders on the far left side of the 
complex s plane. Taking the penalty factor e"^ = lO"* is enough to constrain 
the violation of BNR relations at the level of 1 percent, Further decreasing 
it down to 10"^ does not help much in improving those BNR relations but it 
causes further problems on the stability of the global fit, though the center 
value of the a pole position changes very little. In below we list the fit results 
corresponding to (almost) the minimal x^ we can find: 



Xtot -- 


= Xoo + Xli ± X20 + xInr = 29.7 ± 214.9 ± 41.6 ± 4.4 


Mp-- 


= 757.0 ± OAMeV , Tp = 152.2 ± O.GMeV , 


Mfo 


= 984.5 ± 2.3Mel^ , Tjo = 34.4 ± 6.8MeV , 


M^-- 


= 457 ± l5MeV , T„ = 551 ± 28MeV . 



(21) 



The third pole location in the IJ=00 channel is M3 = 819MeV, Eg = 1846MeV 
with large error bars. The fit results are plotted in Fig. |21 The obtained scat- 
tering lengths and the effective range parameters are listed in table Q where 
we also list several results found in the literature for comparison. From ta- 
ble ^ we find that our results are in general in good agreement with the xPT 
results j2ni and especially the Roy equation analysis [H]. However, since we 
used the flawed low energy data in the IJ=11 channel,''' we will also test in the 
following the manipulated data of 6l from CMD-2. Nevertheless the problem 
of the data in the IJ=11 channel does not seem to propagate into IJ=00 and 20 
channels in our approach, since the constraints from the five BNR relations on 
the IJ=11 channel are rather weak because the magnitude of Tn is suppressed 
in the region governed by the BNR relations. 

In the second approach we use the manipulated data in the IJ=11 channel 



^if wc do not include the so called p' pole the 6} parameter does not agree well with the 
Roy equation analysis. 
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from Ref. j25. The solution is the following: 

xL = Xoo + Xn + Xlo + xInr = 29.8 + 0.02 + 41.6 + 5.9 ; 

Mp = 763.0 ± 0.2MeV , Tp = 139.0 ± 0.5MeV , 

Mp = 984.4 ± 2.7MeV , Tp = 34.6 ± 7.9MeV , 

M„ = 459 ± 3MeV , T^ = 551 ± 23MeV . (22) 

The obtained scattering lengths and the effective range parameters in the 
IJ=00 and 20 channels are found to be in satisfactory agreement with those 
listed in tabled However in the I,J=1,1 channel we find that a\ = 3.54 x 10~^, 
b\ = 4.27 X 10~^, which does not seem to be compatible with the results ob- 
tained from the Wanders sum rule given in Ref. [20] • Nevertheless we find, by 
comparing Eq. ()2H) with Eq. ()22j). that the a pole location changes very little. 

Finally we briefly comment on the use of another group of IJ=20 data by 
Hoogland et al..J22j The quality of solution B of Ref. J22j is not good as it 
gives a very large x^- For solution A, in the combined fit the central value of 
Oq is —4.09 X 10~^. The other outputs are compatible to previous results and 
the fit results in a satisfactory Xbnr ^^^ xlo- O^ ^^^ other side, the mass and 
width of the a resonance, will both increase roughly by 30 MeV. 

A careful analysis was performed to estimate the influence of various theo- 
retical uncertainties to the determination of the a pole, including the penalty 
factor e, various cutoff parameters, different treatment on IJ=11 (and also 
IJ=20) channel data, and whether or not including the third pole. It is found 
that 

M^ = 470 ± bOMeV , T^ = 570 ± bOMeV , (23) 

where the error bars are rather conservative. The results are found in very 
good agreement with the results given in Ref. [B]. Different from Ref. [0], the 
error bars here are mainly caused by the use of different data, rather than by 
analytic continuation. 

To conclude, using the unitarization approach developed by our group and 
combining the constraints from crossing symmetry, it is found that a fit to 
the low energy data gives reasonable results, comparing with Roy equation 
analyses. I2H] The existence of the virtual state pole in the IJ=20 channel 
is reconfirmed. It is also found that crossing symmetry is crucial in correctly 
determining the a pole location. Actually using the present formulae, without 
the constraints from BNR relations the individual fit in the IJ=00 channel 
would give a pole located at Mo- = 542MeV, To- = 546MeV. It is also found 
that the BNR relations provide a stronger correlations between the IJ=00 and 
20 channels and the 11 channel amplitude plays a relatively minor role in the 
BNR relations being used. The fit results end up in general in good agreement 
with the results of Ref. [B]. The estimated value of a pole position as quoted 
by the Ref. j2] is certainly unsatisfactory. It is desirable that newer estimates 
on the pole location begin to converge. 
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